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ABSTRACT 

The  stochastic  forcing  theory  of  Frankignoul  and 
Hasselmann,  1977  is  modified  to  include  a  mixed  layer  model. 
This  enables  the  examination  of  the  interaction  between 
stochastic  heat  flux  or  wind  stress  components  and  the 
annual  period  surface  heat  flux.  The  presence  of  the 
stochastic  heat  flux  component  causes  the  average  sea 
surface  temperature  to  be  0.75  C  higher  then  it  would  be 
with  only  the  annual  period  component.  It  also  delays  the 
time  of  maximum  surface  temperature,  and  it  causes  the 
average  mixed  layer  depth  to  be  ten  meters  shallower.  A 
stochastic  wind  stress  component  applied  to  an  annual  heat 
flux  cycle  produces  a  smaller  sea  surface  temperature  vari- 
ance, but  it  results  in  a  more  rapidly  deepening  and  deeper 
mixed  layer  than  achieved  with  an  annual  cycle  and  constant 
wind  stress.  The  stochastic  forcing  model  shows  that  the 
climatology  of  the  seasonal  thermocline  is  dependent  on 
nonlinear  interactions  between  the  annual  cycle  and 
stochastic  forcing. 
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I.  INTRODUCTION 

The  object  of  this  paper  is  to  show  the  effect  of 
stochastic  surface  heat  flux  forcing  on  the  climatology  of  a 
model  seasonal  thermocline.  It  was  theoretically  demon- 
strated by  Hasselmann  (1976)  that  if  "weather"  forcing  can 
be  considered  to  be  white  noise,  then  the  ocean  acts  as  an 
integrator.  The  result  is  a  red-shifted  response  spectrum. 
Frankignoul  and  Hasselmann  (1977)  developed  a  simple  model 
to  test  this  theory.  Here  we  shall  expand  on  their  result 
by  considering  seasonal  variations  and  mixed  layer  physics. 

Climate  researchers  have  long  been  aware  of  the 
pronounced  variability  associated  with  climatic  records  of 
sea-surface  temperature  (Namias,  1959;  Wyrtki ,  1965).  This 
temporal  variability  spans  a  spectral  range  from  a  low 
frequency  limit  of  the  order  of  ice  ages  (one  cycle  per 
100,000  years)  to  a  high  frequency  limit  on  the  order  of  one 
cycle  per  month. 1  Higher  frequencies  than  one  cycle  per 
month  are  taken  to  be  stochastic  in  this  study,  whereas  mean 
quantities  and  climatologies  will  have  seasonal  or  longer 
frequencies.  An  understanding  of  the  origin  of  this  diverse 
climatic  variability  is  a  primary  goal  of  climate  research. 

Attempts  have  been  made  to  link  climatic  changes  to 
various  external  factors  such  as  solar  activity,  perturba- 
tions in  the  Earth's  orbit,  or  the  increased  atmospheric 
turbidity  associated  with  volcanic  eruptions.  The  primary 
difficulty  with  these  theories  is  that  the  input-response 
relationships  are  not  sufficiently  pronounced  to  be  immedi- 
ately obvious  upon  inspection  of  the  appropriate  time 
series . 


xThe   limits   follow   the  definitions   adopted   in   GARP 
Publication  16,  1975. 
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Recent  interest,  both  commercial  and  tactical,  focuses 
on  the  higher  frequency  sea  surface  temperature  (SST) 
perturbations.  The  1983  El  Nino  (seasonal  frequency  phenom- 
enon) proved  disastrous  to  the  commercial  fishing  industry 
along  the  west  coast  of  the  U.S.  and  South  America. 
Tactically,  the  strategic  defense  of  the  U.S.  is  linked  to 
the  proper  utilization  of  resources  in  the  vicinity  of  SST 
anomalies  and  other  thermal  structures  which  alter  acous- 
tical propagation  through  the  oceans. 

Large  scale  SST-weather  relationships  have  been  discov- 
ered through  the  statistical  analysis  of  extended  time 
series  (Bjerknes,  1966;  Namias,  1969).  Numerical  experi- 
ments have  demonstrated  statistically  significant  modifica- 
tions of  weather  patterns  in  the  tropics  by  critically 
located  SST  anomalies.  These  anomalies  have  been  studied 
through  the  use  of  General  Circulation  Models  (Shukla,  1975; 
Rowntree,  1976)  and/or  Statistical  Dynamical  Models. 

Much  of  the  theory  for  ocean  mixed  layer  modeling 
depends  on  the  validity  of  two  hypotheses.  The  first  is 
that  vertical  mixing  in  the  ocean  mixed  layer  and  entrain- 
ment  at  the  thermocline  occur  as  a  result  of  atmospheric 
forcing:  the  heat  fluxes  across  the  air-sea  surface  and  the 
surface  wind  stress.  The  second  hypothesis  is  that  the 
mechanical  energy  budget  holds  the  key  to  understanding  the 
ocean  mixed  layer  dynamics  (Garwood,  1979). 

The  stochastic  forcing  concept  is  another  approach  to 
understanding  the  ocean  SST  response  to  atmospheric  forcing. 
Lemke  (1977)  applied  the  stochastic  forcing  concept  to  a 
Budyko- Sellers  type  model  to  extend  the  concept  to  a  global 
stochastic  energy  model  applicable  for  time  scales  from  10 
to  10,000  years.  Frankignoul  and  Reynolds  (1983)  examined 
the  statistical  characteristics  of  SST  anomalies  in  the 
presence  of  a  mean  current.  Their  studies  indicate  the 
importance  of  surface  heat  flux  and  wind  mixing  fluctuations 


on  the  production  of  SST   anomalies  from  the  available  North 
Pacific  Experiment  (NORPAX)  data. 

Hasselmann  (1976)  proposed  the  use  of  a  stochastic 
forcing  model  to  study  the  ocean  response  to  atmospheric 
forcing.  He  assumed  that,  because  the  time  scale  of  the 
ocean  response  is  so  much  greater  than  the  time  scale  of  the 
atmospheric  forcing,  the  atmospheric  forcing  can  be  consid- 
ered to  be  represented  as  a  white  noise  time  series.  Via 
this  representation,  he  showed  that  the  slow  climatic  vari- 
ability of  the  ocean  can  be  interpreted  as  an  integral 
response  to  the  random  excitation  provided  by  "weather" 
disturbances . 

In  part  two  of  the  series  on  stochastic  climate  models 
(Frankignoul  and  Hasselmann,  1977),  the  stochastic  forcing 
concept  was  applied  to  a  very  simple  SST  model  utilizing  a 
fixed  depth  mixed  layer,  "copper  plate"  ocean  model.  This 
simplest  model  of  the  ocean  mixed  layer  assumes  a  completely 
mixed,  horizontally  homogeneous  layer  of  depth  h  and  temper- 
ature T  overlying  a  seasonal  thermocline  with  prescribed 
temperature  structure.  If  the  effects  of  salinity  and  hori- 
zontal advection  are  ignored,  the  layer  parameters  T  and  h 
are  governed  by  two  prognostic  equations  of  the  form  (Kraus 
and  Turner,  1967;  Denman,  1973): 


JJp--f,n;T«^,u.R>,  (i-i) 


#--OT-> . 


(1.2) 
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Equation  (1.1)  is  a  heat  balance  equation,  where  f,  denotes 
the  sum  of:  (1)  the  total  flux  of  latent  and  sensible  heat 
through  the  air-sea  interface  (dependent  on  the  air  tempera- 
ture Ta,  humidity  q  and  wind  speed  U  via  standard  bulk  aero- 
dynamic formulas);  (2)  the  short  and  long  wave  radiation  R 
through  the  sea  surface;  and  (3)  the  heat  exchange  with  the 
layer  below.  Equation  (1.2)  is  an  empirical  relation  for 
the  entrainment  of  fluid  from  below  the  mixed  layer  and  is  a 
function  of  wind  speed  and  other  variables. 

Frankignoul  and  Hasselmann  (1977)  simplify  the  model 
further  by  retaining  only  equation  (1.1),  so  that  the  mixed 
layer  depth  remains  constant.  Furthermore,  only  the  contri- 
bution (1)  to  ft  is  included.  Heat  exchange  with  the  atmos- 
phere is  accomplished  by  assuming  that  the  air- sea 
temperature  difference  is  proportional  to  the  north-south 
velocity  V,  via  the  relation  (Ta-T)=KV,  where  K=constant. 
Therefore,  a  wind  from  the  north  will  cool  the  ocean,  and 
conversely  a  wind  from  the  south  will  warm  the  ocean.  Given 
the  above  parameterizations ,  the  rate  of  change  of  SST  is, 


Wl.»^lUl 


df '  V 


The   following  values   were  assigned   to   the  parameters   of 
equation  (1.3): 

K=0.25  (•C/m/s) 

C*  =  10*J 

B=3.0 

/>a  =  1.25  x  10*3  gm  cm"3 

A' =  1.0  gm  cm*3 

c;  =  0.24  cal  gm*Vc)" 

Cj?  =  0.96  cal  gm"'(°C)" 

h=  25  m 
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A  mixed  layer  depth  of  25  m  was  selected  as  representative 
of  lower  mid- latitudes  (30*N)  on  the  basis  of  a  comparison 
of  the  observed  SST  cycle  with  predictions  by  a  "copper 
plate"  model,  (a  value  of  h=100  m  is  used  for  calculations 
at  station  India  59*N,  19°W) .  With  K=0.25,  the  root  mean 
square  (rms)  air-sea  temperature  difference  generated  by  the 
model  was  1.25  C  (Frankignoul  and  Hasselmann,  1977). 

The  atmospheric  model  (Frankignoul  and  Hasselmann,  1977) 
is  a  barotropic  model2  of  quasi-geostrophic  turbulence 
developed  by  Kruse  (1975)  for  numerical  investigations  of 
nonlinear  transfer  within  large  scale  atmospheric  motions  in 
an  equilibrium  state.  It  provides  a  statistically 
stationary  wind  field  from  which  heat  and  momentum  fluxes 
are  inferred  using  bulk  aerodynamic  formulae.  In  this  model 
there  is  no  direct  feedback  from  the  ocean  to  the  atmos- 
phere . 

Kruse' s  model  represents  a  non-divergent,  one  layer 
atmosphere  on  a  /fl  plane  bounded  by  two  latitudes,  and  it  is 
zonally  cyclic.  The  dependent  variable  is  the  stream  func- 
tion which  obeys  the  equation, 


j/j-VV  ♦  JrtC.vY)  V|£  :  forcmg  *yfr  -  r7> 


(1.4) 


where 


^'-^-♦Uo£ 


(1.5) 


2Similar  models  are  discussed  by  Lilly  (1972)  and  Rhines 
(1975). 
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In  equations  (1.4)  and  (1.5),  x  and  y  are  eastward  and 
northward  cartesian  coordinates,  U0  is  a  constant  velocity 
of  5  m/s  and  J  denotes  the  Jacobian  operator  (Frankignoul 
and  Hasselmann,  1977).   The  boundary  conditions  are: 


■« 


-O   Qf  u,o,Lu   ,  d-6) 


and 


Vfx.y.t):  fix  JpL.y.t)    p=l,2,.-         (1.7) 


where,  Lx=Ly=8800  km.  Despite  the  high  degree  of  simplifi- 
cation, the  atmospheric  model  provided  a  statistically 
stationary  wind  field. 

Frankignoul  and  Hasselmann  (1977)  constructed  time 
series  of  the  stochastic  atmospheric  forcing  from  equation 
(1.3)  at  each  gridpoint  and  calculated  the  SST  changes  via 
integration.  To  draw  attention  to  the  evolution  of  low 
frequencies  in  the  SST  fluctuations,  the  SST  time  series  in 
Figure  1.1  have  been  low-pass  filtered.  Despite  the 
simplicity  of  the  coupled  system,  the  white  noise  spectrum 
of  atmospheric  forcing  results  in  a  red-shifted  SST 
response.  Without  a  stabilizing  negative  feedback,  the 
total  SST  variance  is  shown  to  grow  indefinitely  with  time. 
With  a  properly  chosen  negative  feedback,  the  response 
becomes  asymptotically  stationary. 

The  purpose  of  the  present  research  is  to  explore  the 
role  of  seasonality  in  the  ocean  response  to  stochastic 
forcing.    It  is  hypothesized   that  the  seasonal  variability 
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Figure  1.1  Simulated  SST  variations  (without  and  with 
feedback  effect)  and  total  heat  flux  at  x=0,  y=Ly/4. 
The  heat  flux  time  series  is  subsampled  at  5-day 
intervals  and  the  SST  time  series  are  low-passed, 
usine  a  quadratic  Lanczos  filter  (cut-off  frequency  8 
x  10^  Hz  J  (From  Frankignoul  and  Hasselmann,  1977). 


in  the  surface  heat  flux  (Qo),  as  modified  by  a  random 
forcing  component  (Q'),  will  cause  seasonal  variability  in 
the  mean  mixed  layer  depth  (h)  and  in  the  variance  of  mixed 
layer  depth  (h'"*).  Variability  in  h  will  in  turn  modify  the 
mean  temperature  (T)  and  the  temperature  variance  (T'  ). 
This  study  is  part  of  an  investigation  of  increasing 
complexity: 

1)  To  a  constant  wind  stress,    impose  a  seasonal  heating/ 
cooling  (Q0)      cycle  plus  a   stochastic  heat  flux  component 

(Qf); 
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2)  Retain  the  seasonal  heating/cooling  (Q0)  cycle  minus 
the  stochastic  component,  but  add  a  stochastic  component 
to  the  wind  stress;  and 

3)  Consider  the  combined  effect  of  an  annual  heat  flux 
cycle  with  both  a  stochastic  wind  stress  and  heat  flux 
component.  This  step  represents  the  most  realistic  situ- 
ation, but  care  must  be  exercised  to  insure  the  proper 
correlation  between  the  wind  stress  and  surface  heat  flux. 

This  paper  examines  the  ocean  response  to  random  atmospheric 
forcing  as  described  in  investigations  (1)  and  (2)  above. 
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II.  THE  STOCHASTIC  MODEL 

Frankignoul  and  Hasselmann  considered  a  simple  climatic 
system  consisting  of  a  slab  oceanic  mixed  layer  model  which 
was  driven  by  fluxes  of  heat  and  momentum  across  the  air- sea 
interface.  The  governing  equation  (1.3)  represents  a  linear 
approach  to  the  study  of  climate  variability  due  to 
stochastic  forcing.  The  stochastic  model  now  described  is 
not  only  nonlinear  in  its  approach,  but  it  is  designed  to 
simulate  mixed  layer  evolution  on  daily  to  annual  time 
scales  in  response  to  atmospheric  forcing. 

A.   OCEAN  MIXED  LAYER  MODEL  INCLUDING  SEASONAL  THERMOCLINE 

Structure 

A  potentially  significant  improvement  over  previous 
studies  of  stochastic  forcing  lies  in  the  added  complexity 
of  the  ocean  mixed  layer  model.  The  latest  version  of  the 
Garwood  (1977)  one-dimensional  bulk  model  of  the  mixed  layer 
of  the  ocean  was  chosen  for  several  reasons.  First,  the 
model  was  readily  available  and  required  only  minor  modifi- 
cations for  the  task.  Second,  the  model  is  thoroughly 
tested3  and  reproduces  mid- latitude  boundary  layer  dynamics. 
Third,  the  model  produces  solutions  of  the  turbulent  kinetic 
energy  components  and  for  the  terms  of  the  turbulent  kinetic 
energy  budget.  Future  studies  will  more  directly  study  the 
details  of  stochastic  forcing  on  the  turbulence  components. 


3It  is  realized  that  there  are  no  exhaustive  tests  of  a 
model,  but  a  very  successful  study  of  Ocean  Station  Papa 
simulations  performed  in  1981  supports  the  dynamic  consis- 
tency of  the  model  (Garwood  and  Adamec,  1981). 
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The  model  explicitly  parameterizes  a  number  of  mixed- 
layer  processes  ignored  by  other  models.  The  ratio  of  the 
mixed  layer  depth  to  the  Obukhov  length  scale  determines  the 
fraction  of  wind-generated  turbulent  kinetic  energy  avail- 
able for  mixing.  The  viscous  dissipation  is  dependent  on  a 
local  boundary  layer  Rossby  number.  The  model  uses  separate 
equations  for  vertical  and  horizontal  turbulent  kinetic 
energy  (see  Garwood  (1977)  and  Garwood  and  Yun  (1979)  for 
further  details  on  the  mixed  layer  model). 

The  Garwood  model  is  used  here  to  simulate  the  evolution 
of  the  ocean  surface  layer  on  daily  to  annual  time  scales  in 
response  to  atmospheric  forcing.  The  horizontal  gradients 
and  horizontal  advection  are  neglected  in  this  one- 
dimensional  model,  and  the  version  used  in  this  project 
neglects  vertical  diffusion  within  the  thermocline.  In 
addition,  the  salinity  is  held  constant  over  the  water 
column.  The  equation  of  state  is  reduced  to  an  approximate 
proportionality  between  temperature  and  buoyancy. 

In  the  absence  of  a  vertical  advection  velocity,  the  two 
governing  equations  for  the  model  become: 


-fr-^-^ 


(2.1) 


and 


(2.2) 
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where : 

T-mixed  layer  temperature 

Q0 -mean  annual  surface  heat  flux  (positive  downward) 
Q' -random  component  of  the  annual  surface  heat  flux 
AT-temperature  jump  across  the  thermocline 
/Ow-density  of  water 

C,.  -specific  heat  at  constant  pressure  of  water 
h-mixed  layer  depth 
e-entrainment 

Q#  is  prescribed  according  to: 


Q.*A*in(ifi^>     n-i,...,**  ,  (2'3) 


Where  the  amplitude  A  is  equal  to  150  Watts/m1.  Q'  is 
generated  by  means  of  a  uniform  random  number  generator  and 
has  a  maximum  potential  amplitude  equal  to  A.  Therefore, 
the  possible  daily  fluctuation  in  Q*  is  from  plus  to  minus 
A. 

The  mixed  layer  temperature  in  this  model  changes  as  a 
result  of  the  applied  surface  heat  flux  and  entrainment 
mixing  of  colder  water  from  below.  Entrainment  mixing  of 
colder  water  from  below  the  layer  raises  the  potential 
energy  of  the  water  column.  This  energy  for  mixing  is 
supplied  by  the  mixed  layer  turbulence.  Bulk  layer  models 
such  as  this  one  determine  the  rate  of  entrainment  (e)  from 
the  turbulent  kinetic  energy  budget.  In  the  Garwood  model, 
the  time  scale  for  viscous  dissipation  of  turbulent  kinetic 
energy  depends  on  the  intensity  of  the  turbulence,  the  layer 
depth  and  on  the  planetary  rotation.  The  rate  of  entrain- 
ment depends  on  the  vertical  turbulent  kinetic  energy  avail- 
able and   the  temperature  jump  AT   at  the  base  of   the  mixed 
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layer.  Entrainment  becomes  zero  when  the  vertical  turbulent 
kinetic  energy  approaches  zero.  Thus  the  equation  for  the 
rate  of  entrainment  (e)  becomes: 


_  1  


e  s 


m,(ww^(F)  (2.4) 

«q  n  AT 


9 


here 

w'=vertical  turbulent  velocity. 

<E>=mean  mixed  layer  turbulent  kinetic  energy. 

*=thermal  expansion  coefficient. 

g=gravitational  acceleration. 

m  =dimensionless  constant. 

The  mean  mixed  layer  turbulent  kinetic  energy  <!>  is  a  func- 
tion of  both  the  surface  wind  stress  and  the  surface  heat 
flux.  The  relative  distribution  of  the  turbulent  components 
is  particularly  sensitive  to  the  sign  and  magnitude  of  the 
surface  buoyancy  flux.  Positive  buoyancy  flux  produces 
<w'w'>  directly  and  increases  vertical  turbulent  kinetic 
energy  and  entrainment.  Negative  buoyancy  flux  inhibits 
vertical  mixing.  Entrainment  and  mixing  of  water  from  below 
occurs  only  when  there  is  sufficient  turbulent  kinetic 
energy  available  in  the  mixed  layer.  In  this  case  there  is 
active  entrainment  and  mixed  layer  evolution  is  given  by  the 
prognostic  equation  (2.2). 

Under  conditions  of  insufficient  turbulent  kinetic 
energy  for  mixing,  the  mixed  layer  retreats  immediately  to  a 
depth  hr.  This  retreat  depth  is  determined  by  assuming 
steady  state  conditions  for  the  turbulent  kinetic  energy 
budget  and  letting  <w'w'>  approach  zero.  The  retreat  depth 
hr  is  given  by: 
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hr^ff-'jf-^fF) 


(2.5) 


where 

U*  is  the  friction  velocity  at  the  water  surface, 
m ,  and  m3  are  dimensionless  model  constants. 

With  the  prescribed  surface  heat  flux  (Q=Qo+Q' )  and  the 
friction  velocity  U-,  the  model  calculates  <E>  and  <w' w' > 
and  the  mixed  layer  temperature  and  depth  from  (2.1)  to 
(2.5). 

1.   Reduction  to  "Copper  Plate"  Form 

The  Garwood  (1977)  model  is  easily  converted  to  a 
simple  slab  or  "copper  plate"  model  as  used  by  Frankignoul 
and  Hasselmann  (1977).  This  is  accomplished  by  setting  the 
mixed  layer  depth  (h)  equal  to  some  constant  value.  The 
result  is  that  equation  (2.2)  becomes: 


#■ 


(2.6) 


Specifically,    this  implies   infinite  thermal   conductivity 

(strong  mixing)  between  Z=-h  and  Z  =  0,  but  zero  conductivity 

(no  mixing)    below  this.    Given  this  condition,   equation 
(2.1)  reduces  to: 


iF-i^fr  • 


Equation  (2.7)    is  effectively  the   form  of   equation  (1.3) 
that  Frankignoul  and  Hasselmann  used. 
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The  primary  difference  between  (2.7)  and  (1.3)  lies 
in  the  determination  of  the  heat  transfer  between  the  ocean 
and  atmosphere.  Frankignoul  and  Hasselmann  (1977)  parame- 
terize the  heat  transfer  as  a  function  of  the  wind  direc- 
tion, while  Garwood  (1977)  specifies  the  surface  heat  flux 
directly.  However,  both  (2.7)  and  (1.3)  will  produce  equiv- 
alent stochastic  forcing.  The  advantage  of  (2.7)  is  that 
the  rate  of  change  of  SST  is  now  dependent  upon  both  the 
annual  cycle  and  the  random  forcing,  and  facilitates  later 
application  with  a  predictive  equation  for  mixed  layer 
depth. 

Reducing  Garwood's  model  to  a  simple  slab  model 
permitted  the  recreation  of  the  results  of  Frankignoul  and 
Hasselmann  (1977)  for  the  case  without  negative  feedback. 
These  results  will  be  discussed  further  in  Chapter  III. 

B.   ATMOSPHERIC  FORCING  INCLUDING  AN  ANNUAL  CYCLE 

The  ocean  model  is  driven  by  an  annual  period  sinusoidal 
surface  heat  flux  (QD),  a  stochastically  varying  heat  flux 
(Q' )  and  a  constant  wind  stress.  For  computational  conven- 
ience a  model  year  is  taken  to  be  360  days.  The  model  is 
initialized  at  "spring  transition"  corresponding  approxi- 
mately to  the  vernal  equinox  with  a  mixed  layer  depth  of  100 
meters,  a  mixed  layer  temperature  of  6  C  and  an  underlying 
homogeneous  water  column  of  5  C.  The  magnitude  of  the 
annual  period  surface  heat  flux  was  chosen  so  that  the 
surface  temperature  would  rise  approximately  8  C  and  the 
seasonal  thermocline  would  shallow  to  approximately  50 
meters  at  the  times  of  maximum  SST  and  maximum  downward 
surface  heat  flux  respectively.  The  necessary  amplitude  of 
the  annual  heat  flux  to  achieve  this  was  150  Watts/miin  the 
case  of  a  1  dyne/cm*  wind  stress.  The  salinity  is  specified 
as  a  constant  33  %©  and  it   has  no  dynamical  significance  in 
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this  study.  The  initialization  described  approximates 
conditions  at  Ocean  Station  Papa  at  50°N,  145°W.  To  allow 
the  model  to  reach  a  statistical  equilibrium,  the  model  is 
actually  integrated  to  720  days,  and  the  second  model  year 
is  used  for  calculations  and  graphic  display. 

The  stochastic  forcing  is  supplied  by  a  random  number 
generator.  A  known  but  different  seed  value  is  specified 
for  each  case,  which  allows  duplication  of  the  random 
forcing  series.  The  stochastic  forcing  is  applied  to  the 
surface  heat  flux  such  that  the  random  forcing  can  on  any 
day  equal  (positively  or  negatively)  the  amplitude  of  the 
annual  surface  heat  flux  cycle. 

C.   CALCULATION  OF  STATISTICAL  VARIABLES 

At  this  point  it  becomes  necessary  to  define  the  nomenc- 
lature to  be  used  in  discussions  of  the  model  results. 
Model  outputs  are  sea-surface  temperature,  mixed  layer 
depth,  the  average  temperature  variance  T'*  and  average 
mixed  layer  depth  variance  h'*.  The  variables  T0  and  h0 
represent  the  ocean  response  in  temperature  and  mixed  layer 
depth  to  the  annual  period  sinusoidal  heat  flux  cycle 
without  the  stochastic  component.  The  differences  T-T0  and 
h-h0,  between  the  mean  surface  temperature  or  mixed  layer 
depth  with  random  forcing  and  the  mean  surface  temperature 
or  mixed  layer  depth  without  random  forcing  provide  insight 
into  the  effects  of  the  stochastic  forcing  theory. 

An  overbar  denotes  an  ensemble  mean  of  the  daily  values 
of  the  variable  divided  by  the  number  of  cases  summed. 
Thus,  the  output  statistics  represent  daily  averages  over 
some  climatological  time  frame.  The  variance  of  temperature 
and  mixed  layer  depth  can  be  computed  from  the  following 
equations  obtained  after  performing  a  Reynolds 
decomposition. 
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t":CtV(T)2    ,  (2.8) 


and 


'i    "I.      -.1  (2-9) 


h   :(h')  -(h)' 


Thus,  the  variances  are  the  result  of  the  difference  between 
mean  quantities. 
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III.  THE  STOCHASTIC  MODEL  EXPERIMENTS 

A.   PROCEDURES 

1 .  General  Method 

The  stochastic  forcing  is  applied  to  101  cases, 
which  corresponds  to  101  model  years.  The  first  case  or 
model  year  represents  the  nonstochastic  case  that  determines 
the  daily  value  of  T0  and  h0 .  The  next  100  cases  contain  a 
daily  random  forcing  component  to  the  annual  surface  heat 
flux  cycle.  From  these  cases,  the  averages  and  variances  of 
the  temperatures  and  mixed  layer  depths  are  calculated.  The 
average  temperature  and  mixed  layer  depth  are  computed  daily 
via  a  simple  summation  over  the  100  stochastically  forced 
cases.  For  the  variance  computation  see  equations  (2.8)  and 
(2.9). 

2.  Reproducing  Frankignoul  and  Hasselmann' s  Result 


Frankignoul  and  Hasselmann  (1977)  applied  the 
stochastic  forcing  concept  to  a  simple  slab  or  "copper 
plate"  ocean  model  as  given  by  equation  (1.3).  The  results 
obtained  from  the  Garwood  model  as  in  equation  (2.7)  confirm 
the  hypothesis  of  Hasselmann  (1976)  that  the  ocean  acts  as 
an  integrator  of  the  white  noise  atmospheric  forcing  and 
produces  a  red-shifted  response  spectrum,  that  is  the 
spectral  response  of  the  model  is  shifted  toward  the  lower 
frequencies,  and  appears  geophysical  in  nature  (Figure  1.1). 
Frankignoul  and  Hasselmann  further  showed  that,  without  a 
stabilizing  negative  feedback,  the  total  SST  variance  grows 
indefinitely  with  time,  as  suggested  by  Hasselmann  (1976). 

Figure  3.1  displays  the  daily  model  response  to  the 
sinusoidal   annual   heat   flux   cycle   without   stochastic 
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forcing.  It  can  be  seen  that  the  annual  surface  heat  flux 
cycle  begins  at  vernal  equinox,  and  adds  heat  to  the  mixed 
layer  with  the  maximum  heating  rate  occurring  at  90  days. 
Heat  continues  to  be  added  to  the  mixed  layer  for  an  addi- 
tional 90  days,  and  at  day  180  the  cycle  reverses  and  heat 
is  withdrawn  from  the  ocean,  with  the  maximum  withdrawal 
rate  on  day  270.  Heat  continues  to  be  removed  from  the 
ocean  until  the  cycle  is  complete  on  day  360. 

Although  the  one-dimensional  model  is  designed  so 
that  the  mixed  layer  temperature  will  rise  approximately  8  C 
with  realistic  mixed  layer  depths,  it  is  evident  that  fixing 
the  mixed  layer  depth  reduces  the  ocean  temperature  response 
to  only  approximately  3.5  C.  Since  the  mixed  layer  depth  is 
constant  there  is  no  interaction  between  the  minimum  depth 
at  the  time  of  maximum  surface  heat  flux.  Instead,  the  heat 
input  is  spread  over  the  entire  100  m,  and  the  temperature 
increase  is  only  about  half  of  the  expected  value. 

The  effects  of  the  stochastic  forcing  for  a  typical 
case  can  be  seen  in  Figure  3.2.  The  SST  would  statistically 
be  expected  to  be  larger  than  T6  one  half  of  the  time  and 
smaller  than  T0  the  other  half.  In  this  case  the  SST  is 
slightly  smaller  than  T„. 

Figure  3.3  shows  the  significant  temperature  statis- 
tics calculated  from  the  model.  Hasselmann  (1976)  and 
Frankignoul  and  Hasselmann  (1977)  showed  that  the  SST  vari- 
ance in  the  absence  of  negative  feedback  will  grow  indefi- 
nitely. The  upper  graph  of  Figure  3.3  confirms  this  result. 
The  variance  of  the  temperature  perturbation  can  clearly  be 
seen  to  exhibit  a  ramp  type  growth  over  the  period  of  the 
experiment.  Also,  there  is  no  statistically  significant 
difference  between  the  average  temperature  with  random 
forcing  and  the  average  temperature  response  to  the  annual 
cycle  alone,  T-T„.  If  the  experiment  was  repeated  an 
infinite  number  of  times,   the  values  of  T-T„  would  be  zero. 
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Figure  3.2  Typical  year  with  stochastic  forcing.  a) 
The  same  annual  heat  flux  cycle  shown  in  Figure  3.1, 
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b)  The  SST  over  the  model  year.  The  mixed  layer  depth 
h=100  m. 
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Figure  3.3  Statistical  temperature  variables.  a)  The 
mean  temperature  variance  of  100  stochastically  forced 
cases.   fa)  Plot  of  T-Te . 
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This  result  only  holds  for  the  case  where  the  mixed  layer 
depth  is  fixed,  which  is  a  physically  unrealistic  situation. 
From  the  above  results,  it  can  be  said  that  the 
modified  Garwood  model  in  slab  form  correctly  duplicates  the 
linear  time  rate  of  change  of  SST  developed  by  Frankignoul 
and  Hasselmann  (1977).  Because  the  Garwood  model  is  not 
linear  and  the  forcing  is  completely  specified,  it  is  not 
possible  to  compute  an  adequate  negative  feedback  term, 
without  first  determining  the  climatology,  to  confirm 
further  studies  made  by  Frankignoul  and  Hasselmann.  In  the 
next  section,  we  explore  the  role  of  seasonality  in  the 
ocean  response  to  stochastic  forcing  using  the  one- 
dimensional  model  with  entrainment  into  the  mixed  layer. 

3 .   Stochastic  Forcing  with  Mixed  Layer  Physics 

For  this  experiment,  the  governing  equations  are 
(2.1),  (2.2),  and  the  rate  of  entrainment  is  given  by  equa- 
tion (2.4).  This  study  presents  a  significant  departure 
from  the  slab  models  used  by  Hasselmann  (1976)  and 
Frankignoul  and  Hasselmann  (1977),  because  the  mixed  layer 
temperature  is  no  longer  independent  of  the  mixed  layer 
depth.  Given  a  constant  wind  stress,  the  mixed  layer  depth 
will  shallow  during  periods  of  strong  downward  surface  heat 
flux,  which  also  causes  an  increase  in  mixed  layer  tempera- 
ture. Thus  on  an  annual  basis,  the  effects  of  the  heat  flux 
and  shallow  mixed  layer  depth  will  reinforce  each  other  and 
generate  larger  perturbations  in  temperature.  The  atmos- 
pheric forcing  and  model  initializations  will  remain  as 
stated  for  the  previous  study. 

Figure  3.4  displays  the  model  response  to  the  annual 
period  atmospheric  surface  heat  flux  without  a  stochastic 
forcing  component.  As  stated  earlier  the  atmospheric 
forcing  amplitude  is  designed  to  allow  the  SST  to  increase 
by  approximately  8  C  and  the  mixed  layer  depth  to  shallow  to 
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Figure  3.4  Annual  cycle  without  stochastic  forcing, 
a)  The  daily  SST  response  to  the  annual  cycle.  b)  the 
daily  mixed  layer  depth  response  to  the  annual  cycle. 
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approximately   50  m.     These   conditions   are  met   in   this 
version  of  the  model. 

The  SST  reaches  a  maximum  at  approximately  day  170 
which  indicates  a  80  day  lag  from  the  point  of  maximum  down- 
ward heat  flux.  As  the  second  model  year  of  each  case  is 
used  for  calculations,  the  mixed  layer  depth  (Figure  3.4)  at 
vernal  equinox  is  slightly  shallower  than  the  100  m  used  to 
initialize  the  model.  The  depth  continues  to  shallow  and 
reaches  minimum  values  at  about  day  100.  For  the  remainder 
of  the  model  year,  the  mixed  layer  depth  gradually  deepens 
until  approximately  day  350.  From  day  310  to  day  350,  the 
mixed  layer  depth  has  achieved  the  greatest  depth  that 
turbulent  mixing  will  allow,  and  it  is  entraining  the  top  of 
the  permanent  thermocline.  During  the  last  10  days,  the 
mixed  layer  begins  to  shallow  toward  the  initial  value. 

Upon  comparing  Figures  3.4  and  3.1,  it  is  apparent 
that  the  maximum  SST  is  greater  when  the  mixed  layer  depth 
is  allowed  to  vary.  As  the  atmosphere  adds  heat  to  the 
ocean,  the  mixed  layer  becomes  more  stable,  and  the  poten- 
tial energy  of  the  system  increases.  Increasing  the  poten- 
tial energy  of  the  mixed  layer  is  at  the  expense  of  the 
turbulent  kinetic  energy,  which  causes  the  mixed  layer  depth 
to  retreat  according  to  equation  (2.5).  The  initial  shal- 
lowing of  the  mixed  layer  corresponds  to  the  period  of 
reduced  surface  cooling  in  late  winter  and  early  spring  when 
entrainment  ceases  because  of  the  reduced  buoyancy  produc- 
tion of  turbulent  kinetic  energy  and  a  relative  increase  in 
the  dissipation  of  turbulent  kinetic  energy.  Equation  (2.1) 
can  be  reduced,  for  this  case  with  forced  heat  fluxes,  to  a 
ratio  of  the  atmospheric  heat  flux  and  the  mixed  layer 
depth.  As  the  time  rate  of  change  in  temperature  is 
inversely  proportional  to  the  mixed  layer  depth,  the  shal- 
lowing of  the  mixed  layer  coupled  with  the  downward  heat 
flux  produces  a  rapid  increase  in  temperature  with  time. 
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Figure  3.5  is  for  a  typical  model  year  under  the 
influence  of  both  seasonal  and  stochastic  forcing.  The 
daily  random  heat  flux  component  Q'  is  identical  to  the 
forcing  in  Figure  3.2.  In  all  of  the  graphs  the  annual 
cycle  is  evident.  The  first  significant  result  can  be  seen 
in  the  graph  of  mixed  layer  depth  for  a  typical  year  in 
Figure  3.5.  It  is  evident  that  the  daily  variations  in 
mixed  layer  depth  can  be  as  great  as  50  m  during  portions  of 
the  year.  Therefore,  a  linear  model  of  stochastic  forcing 
with  fixed  layer  depth  is  unrealistic.  The  temperature 
perturbation  is  dependent  on  the  interaction  between  the 
annual  plus  stochastically  forced  heat  fluxes  and  the 
perturbations  in  the  mixed  layer  depth. 

Figure  3.6  shows  the  temperature  statistics  calcu- 
lated from  the  100  cases  in  the  experiment.  There  is  a 
phase  shift  associated  with  the  time  of  maximum  SST,  with 
the  maximum  SST  occurring  about  two  weeks  later  for  the 
stochastically  forced  cases.  The  difference  between  this 
result  and  that  obtained  for  the  slab  model,  can  be  attrib- 
uted to  the  combined  effects  of  a  changing  mixed  layer  depth 
and  the  evolving  SST.  This  curve  does  retain  a  strong 
annual  cycle. 

The  temperature  variance  curve  of  Figure  3.6  retains 
the  increased  amplitude  throughout  the  year  observed  in 
Figure  3.3  for  the  slab  model,  but  the  influence  of  the 
annual  cycle  is  clearly  evident  and  produces  a  1  C  maximum 
temperature  variance.  The  growth  visible  in  the  variance 
graph  is  a  result  of  the  lack  of  any  negative  feedback 
applied  to  the  model. 

To  further  delineate  the  temperature  differences, 
the  graph  of  T-T«  reveals  that  the  stochastically  forced 
average  temperature  is  approximately  0.75  C  higher  than  the 
temperature  derived  from  the  annual  cycle  alone.  This 
result  is   significant.    Both  curves   have  a   strong  annual 
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cycle,   which  indicates  that   the  seasonal  variations  cannot 
be  ignored  when  studying  SST  perturbations. 

Figure  3.7  shows  the  mixed  layer  depth  statistics 
analogous  to  those  computed  for  temperature.  It  is  apparent 
that  the  average  mixed  layer  depth  with  stochastic  forcing 
is  smaller  than  the  mixed  layer  depth  produced  by  the  annual 
cycle  alone.  A  considerable  phase  shift  in  the  minimum 
mixed  layer  depth  occurs.  The  stochastically  forced  depth 
minimum  occurs  approximately  one  month  later  than  the 
minimum  in  h0 .  It  is  this  phase  difference  which  produces  a 
corresponding  phase  shift  in  the  temperature  statistic, 
T-T0.  As  a  result  of  this  phase  shift  a  second  harmonic 
becomes  evident  in  the  h-h„  curve.  The  final  effect  of  note 
pertains  to  the  mixed  layer  depth  variance.  The  mixed  layer 
exhibits  greater  fluctuations  during  the  shallowing  regime 
than  the  deepening  regime.  This  phenomenon  can  be  explained 
through  an  examination  of  the  turbulent  kinetic  energy 
budget.  In  the  spring  the  near-zero  heat  flux  Q„,  allows 
the  mixed  layer  depth  to  be  controlled  by  Q' ,  the  random 
heat  flux  component.  At  this  time,  h'  will  show  a  large 
variation  in  proportion  to  Q' .  In  addition,  there  is  a 
small  temperature  jump  (AT)  across  the  seasonal  thermocline 
and,  as  can  be  seen  in  equations  (2.2)  and  (2.4),  the  time 
rate  of  change  of  mixed  layer  depth  is  inversely  propor- 
tional to  AT.  In  the  fall,  the  temperature  jump  across  the 
seasonal  thermocline  is  very  large.  This  greatly  restricts 
mixed  layer  depth  fluctuations.  Thus  the  shallowing  regime 
has  a  greater  mixed  layer  depth  variance  than  does  the  deep- 
ening regime.  The  curve  of  h-h0  is  always  negative,  except 
for  a  random  phenomenon  occurring  around  day  20,  which  indi- 
cates that  the  stochastically  forced  mixed  layer  is  on 
average  20  m  shallower  than  the  mixed  layer  produced  solely 
by  the  annual  cycle. 


35 


RVERRGE  H  PRIME  SQURRED 


0.0   30.0    60.0    90.0    120.0  150.0  180. 0  210.0  240.0  270.0  300.0  330.0  360. 0 


TIME   (DAYS) 


HBRR  -  HZERO 


u 

U  o 


0-  o 

U    ri 

°8 


NO.   DF  YERRS 

100) 

i           i 

i  ^ 

1 

r       T 

i 

i 

i 

0.0   30.0    60. 0    90.0    120.0  150.0  180.0  210.0  240.0  270.0  300.0  330.0  360.0 

TIME   (DAYS) 


Figure  3.7   Mixed. layer  depth   statistics-full  model. 


a)   The   variance  in   the  mixed   layer  depth   over  the 

The  difference  between  the  average 
when  randomly  forced  and  the  mixed 
is  the  product  of   the  annual  cycle 


experiment.  m   b) 


mixed  layer  depch 
layer  depth  which 
alone,  h-ho. 


36 


B.   INTERPRETATIONS 

1.   "Slab  Model"  Experiment 

With  a  constant  depth  mixed  layer  the  results 
obtained  confirm  those  of  previous  researchers.  The  ocean 
does  respond  to  white  noise  atmospheric  forcing  with  a 
"red-shifted"  spectrum.  The  annual  cycle  affects  only  the 
average  temperature,  not  the  variance  in  temperature. 
However,  the  temperature  variance  does  display  the  indefi- 
nite growth  predicted  by  Hasselmann,  since  there  is  no 
negative  feedback. 

The  assumptions  used  in  this  experiment  produce  a 
linear  equation  for  the  change  in  T'  .  There  are  no 
temperature-mixed  layer  depth  interactions.  This  simple 
model  eliminates  the  possible  contribution  of  the  applied 
annual  period  heat  flux  forcing.  The  results  would  be  the 
same  if  the  random  surface  heat  flux  component  were  applied 
to  a  constant  surface  heat  flux.  Therefore,  the  initial 
model  assumptions  of  Hasselmann  and  Frankignoul  and 
Hasselmann  oversimplify  the  problem.  To  produce  more  real- 
istic results  requires  the  added  complexity  associated  with 
the  nonlinear  dependence  between  mixed  layer  depth  and  SST. 

For  simplicity,  assume  that  the  time  rate  of  change 
of  SST  is  given  as  follows: 


41  .  JL  (3-1) 


By  means  of  a  Reynolds  decomposition  it  can  be  shown  that  it 
is  possible  to  linearize  equation  (3.1)  only  if  I  h'  I  <<  h. 
Expanding  equation  (3.1)  to  include  the  average  and  pertur- 
bation terms  results  in  the  following  expression 


37 


AT  .  AT '  .    Q  *Q     -  Q    j  Q_ 


A  binomial  expansion  gives 


h         h     h  h 


Therefore  the  approximate  form  of  equation  (3.1)  becomes 


it  *  At    "      u      l        u 


If  Ih'l  <<  h,  equation  (3.4)  reduces  to 


_AI  f  AT.  .  _Q_  ,  _qL      . 


Averaging  each  term  produces 


(3.4) 


(3.5) 


-ff  ♦  -M-  r  -fl-  ♦  -S-    .  (3.6) 

o 

Subtracting  (3.6)  from  (3.5)  results  in  the  following  linear 
equation  for  the  change  in  T' 

This   equation  is   the  functional   form  used  by  Hasselmann 
(1976).    It   is  evident   that  the   temperature  perturbation 
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4t   '    L     •  (3.7) 


will  be  the  integral  of  the  right  hand  side.  Thus  the  ocean 
acts  as  an  integrator  of  the  random  forcing.  It  is  equally 
evident  that  the  annual  heat  flux  cycle  plays  no  part  in  the 
determination  of  the  temperature  perturbation. 

2.   One -Dimensional  Model  Experiment 

Suppose  Ih'l   is  not  much   less  then  h.    Expanding 
equation  (3.4)  gives 


XL. ,  XL,  =  _Q_  .  _aC  .  JL  .  .qV  (3.8) 


Averaging  each  term  produces 


(3.9) 


xl  ♦  .5' ..  _a. .  M-  ♦  -au  -  4^ 


Subtracting  (3.9)  from  (3.8)  gives 


/  » 


XL-  _Q_  .  J2_h_  .  _Q_b_  .  _Qji' 


(3.10) 


Here  the  time  rate  of  change  of  T'  depends  on  the  correla- 
tion between  the  random  forcing  component  and  the  perturba- 
tion mixed  layer  depth  and  on  the  relationship  between  the 
average  heat  flux  and  the  mixed  layer  depth  perturbation. 
If,  as  in  the  second  experiment,  Ih'l  is  not  much  smaller 
than  h,  a  linear  relationship  for  T'  is  not  possible. 
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The  results  of  these  two  experiments  clearly  demon- 
strate the  complex  relationships  existing  in  a  single  phase 
of  ocean  response.  The  linear  approach  to  stochastic 
forcing  does  bring  out  some  basic  responses,  but  ignores 
some  very  significant  interactions.  The  ocean  does  respond 
as  an  integrator  of  the  random  forcing  and  the  variance  in 
SST  does  grow  indefinitely  when  no  feedback  is  applied,  but 
the  strong  influence  of  the  annual  cycle  is  missed  when  a 
slab  model  is  used.  In  addition,  the  dynamics  of  the  slab 
model  are  such  that  the  average  of  the  stochastically  forced 
climatology  is  equivalent  to  the  annual  cycle  climatology. 
This  is  not  realistic.  The  two  climatologies  are  not  iden- 
tical. The  random  forcing  component  imparts  enough  extra 
energy  to  the  ocean  mixed  layer  to  raise  the  temperature  by 
0.75  C  during  the  summer,  and  it  delays  the  temperature 
maximum  and  depth  minimum  until  later  in  the  year.  These 
effects  are  important. 

C.   STOCHASTIC  WIND  STRESS  EXPERIMENT 

The  next  step  in  the  investigation  of  the  stochastic 
forcing  theory  is  to  examine  the  model  response  to  a  wind 
stress  composed  of  a  mean  component  and  a  random  perturba- 
tion. Previous  investigations  into  the  importance  of  the 
wind  stress  in  the  development  of  SST  anomalies  indicate 
this  is  an  important  forcing  mechanism  (Camp  and  Elsberry, 
1978;  Frankignoul  and  Reynolds,  1983). 

1.   General  Method 

A  modified  version  of  the  Garwood  (1977)  model 
described  in  the  previous  sections  is  used  for  this  investi- 
gation. The  model  is  modified  such  that  there  is  no  longer 
a  stochastic  heat  flux  component  Q' ,  but  the  annual  period 
heat  flux  cycle  remains.    The   kinematic  wind  stress  U*   is 
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decomposed  into  a  mean  component  U*  and  a  stochastic  compo- 
nent (U*2)'.  The  mean  kinematic  wind  stress  is  1.0  cm*/s* 
and  the  random  component  is  allowed  a  daily  variation  of 
plus  or  minus  0.9  cmz/sz.  *n  all  other  respects  the  proce- 
dure is  the  same  as  that  used  in  the  stochastic  heat  flux 
investigations . 

2.   Stochastic  Wind  Stress  with  Seasonal  Thermocline 

This  experiment  isolates  the  model  response  to 
stochastically  forced  wind  stress  and  a  simple  sinusoidal 
annual  surface  heat  flux  while  retaining  a  variable  mixed 
layer  depth.  As  in  previous  investigations,  the  first  case 
possesses  no  stochastic  forcing  to  determine  the  base  or 
zero  state.  The  model  is  then  run  100  years  and  the  output 
averaged  to  approximate  climatology.  The  temperature  vari- 
ance (Figure  3.8)  does  not  appear  to  display  an  indefinite 
growth  in  spite  of  the  fact  that  there  is  no  feedback 
applied  to  the  model.  Such  growth,  although  still  present, 
is  insignificant  at  the  scale  of  the  figure.  Additionally, 
the  maximum  temperature  variance  occurs  about  day  112, 
approximately  40  days  prior  to  the  maximum  average  SST.  The 
magnitude  of  the  temperature  variance  is  much  smaller  than 
for  the  stochastic  heat  flux  experiment.  The  annual  cycle 
is  evident  in  both  curves  of  Figure  3.8.  The  mean  tempera- 
ture response  to  stochastic  wind  forcing  is  smaller  than  the 
temperature  response  to  the  annual  cycle  alone.  It  can  be 
seen  that  the  temperature  of  the  stochastically  forced  cases 
is  about  0.25  C  smaller  than  Te ,  the  temperature  response  to 
the  annual  cycle  alone. 

Figure  3.9  shows  the  mixed  layer  depth  statistics 
for  the  random  wind  stress  experiment.  The  magnitude  of  the 
variance  of  the  mixed  layer  depth  in  the  shallowing  regime 
is  greater  than  it  was  for  the  heat  flux  experiment  (Figure 
3.7).    During   the  deepening  regime,    the  variance   in  h' 
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approaches  zero.  This  phenomenon  can  be  explained  through 
an  examination  of  the  turbulent  kinetic  energy  budget  (see 
discussion  for  Figure  3.7).  Examination  of  the  lower  graph 
of  Figure  3.9  reveals  that  the  stochastically  forced  mixed 
layer  depth  is  about  15  m  smaller  than  h  during  the  shal- 
lowing regime.  In  comparison  with  Figure  3.7,  it  is  evident 
that  the  mean  mixed  layer  depth  with  stochastic  wind  stress 
forcing  is  about  5  m  smaller  than  with  stochastic  heat  flux 
forcing.  During  the  deepening  regime,  however,  the  mean 
mixed  layer  with  stochastic  wind  forcing  is  actually  deeper 
than  h  the  model  response  to  the  annual  cycle.  Therefore, 
the  stochastic  wind  stress  component  causes  the  mixed  layer 
to  deepen  more  rapidly  through  additional  generation  of 
turbulent  kinetic  energy  that  is  not  completely  compensated 
by  dissipation. 

In  conclusion,  it  can  be  stated  that  the  stochastic 
forcing  concept  when  applied  to  the  wind  stress  U**  does 
affect  the  SST  response  of  the  ocean  model  by  as  much  as 
0.4  C.  The  stochastic  wind  stress  forcing  reduces  SST  on 
the  whole  because  of  enhanced  mixing.  The  random  wind 
stress  forcing  becomes  an  important  source  of  turbulent 
kinetic  energy  during  the  deepening  regime  in  the  late 
summer  and  early  winter.  The  mixed  layer  deepens  more 
rapidly  and  achieves  a  greater  depth  than  without  the 
stochastic  wind  stress  component. 
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IV.  SUMMARY  AND  RECOMMENDATIONS 

The  Hasselmann  stochastic  forcing  theory  was  verified 
for  a  slab  model,  and  then  it  was  modified  to  include  a 
mixed  layer  depth  model.  This  enabled  the  exploration  of 
the  interaction  between  stochastic  and  annual  period  forcing 
of  the  ocean.  The  simple  linear  prediction  equation  for  T' 
that  was  developed  by  Frankignoul  and  Hasselmann  was  no 
longer  applicable  because  of  nonlinear  interactions  between 
h'  and  T'.  The  mixed  layer  depth  perturbations  are  signifi- 
cant in  comparison  with  the  average  mixed  layer  depth.  The 
temperature  perturbation  is  a  function  not  only  of  the 
stochastic  heat  flux  component,  but  also  of  the  interaction 
between  the  annual  heat  flux  cycle  Qc  ,  the  stochastic 
forcing  component  Q' ,  and  the  mixed  layer  depth  perturba- 
tion. Of  considerable  significance  is  the  fact  that  with 
the  stochastic  forcing  the  climatological  SST  response  is 
0.75  C  higher  than  the  corresponding  ocean  response  to  the 
annual  cycle  alone.  This  is  consistent  with  the  climatolo- 
gical mixed  layer  depth  being  smaller  when  random  forcing  is 
applied.  Furthermore,  the  minimum  mixed  layer  depth  occurs 
almost  one  month  later  in  the  year  than  the  ocean  response 
to  the  annual  cycle  alone.  Therefore,  the  stochastic 
forcing  component  of  the  atmospheric  heat  flux  greatly 
influences  the  climatology  of  both  the  SST  and  mixed  layer 
depth. 

In  confirmation  of  Elsberry  and  Camp  (1978),  we  show 
that  the  wind  stress  with  a  stochastic  component  produces  a 
cooler  (by  0.25  C)  ocean  surface  than  would  be  observed 
without  the  random  component.  This  is  because  the  mixed 
layer  is  deeper  on  the  average  due  to  the  generation  of 
additional  turbulent  kinetic  energy.  Interestingly,  however, 
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the  climatological   mixed  layer   depth  is   reduced  in  value 
during  much  of  the  cooling  season. 

Recommended  future  effort  will  be  to  apply  the 
stochastic  forcing  simultaneously  to  the  wind  stress  and  the 
heat  flux.  This  represents  the  most  realistic  situation, 
but  care  must  be  exercised  to  insure  the  proper  correlation 
between  the  random  forcing  of  each  variable,  because  wind 
stress  and  heat  flux  are  not  independent.  This  is  the  first 
step  to  analyzing  the  effect  of  stochastic  forcing  on  the 
individual  turbulence  components. 

The  SST  response  of  the  ocean  on  climatological  time 
scales  is  dependent  on  both  annual  cycles  and  stochastic 
forcing  components.  The  interactions  are  complex  and 
nonlinear.  Only  by  examining  the  parts  is  it  possible  to 
understand  the  whole,  but  care  must  be  exercised  to  insure 
that  the  simplifications  used  do  not  misrepresent  the  actual 
interactions . 
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